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Abstract 



A general relationship is presented between the statistics of thermal radia- 
tion from a random medium and its scattering matrix S. Familiar results for 
black-body radiation are recovered in the limit S — * 0. The mean photocount 
n is proportional to the trace of 1 — S ■ , in accordance with Kirchhoff's 
law relating emissivity and absorptivity. Higher moments of the photocount 
distribution are related to traces of powers of 1 — S ■ <ív, a generalization of 
Kirchhoff's law. The theory can be applied to a random amplifying medium 
(or "random làser" ) below the làser threshold, by evaluating the Bose-Einstein 
function at a negative temperature. Anomalously large fluctuations are pre- 
dicted in the photocount upon approaching the làser threshold, as a con- 
sequence of overlapping cavity modes with a broad distribution of spectral 
widths. 

To appear in: Diffuse Waves in Complex Media, edited by J.P. Fouque, 
NATO ASI Series (Kluwer, Dordrecht, 1999). 
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I. INTRODUCTION 



The name "random làser" made its appearance a few years ago M, in connection with 
experiments on amplifying random media 0. The concept goes back to Letokhov's 1967 
proposal to use a mirrorless làser as an optical frequency Standard || . Làser action requires 
gain and feedback. In any làser, gain results from stimulated emission of radiation by àtoms 
in a non-equilibrium state. The random làser differs from a conventional làser in that the 
feedback is provided by múltiple scattering from disorder rather than by confinement from 
mirrors. Because of the randomness, there is no geometry-dependent shift of the làser line 
with respect to the atòmic transition frequency (hence the potential as a frequency Standard). 
Stellar atmospheres may form a naturally occuring realization of a random làser |4| . 

Possible applications as "paint-on làsers" [BJ have sparked an intensive experimental and 
theoretical investigation of the interplay of múltiple scattering and stimulated emission. The 
tòpic has been reviewed by Wiersma and Lagendijk (see Ref. || and these Proceedings). 
A particularly instructive experiment J7| was the demonstration of the narrowing of the 
coherent backscattering cone as a result of stimulated emission below the làser threshold. 
This experiment can be explained within the framework of classical wave òptics. Wave 
òptics, as opposed to ray òptics, because coherent backscattering is an interference effect. 
Classical òptics, as opposed to quantum òptics, because stimulated emission can be described 
by a classical wave equation. (What is needed is a dielectric constant with a negative 
imaginary part.) 

In a recent work || we went beyond classical òptics by studying the photodetection 
statistics of amplified spontaneous emission from a random medium. Spontaneous emission, 
as opposed to stimulated emission, is a quantum optical phenomenon that can not be de- 
scribed by a classical wave equation. In this contribution we review our theory, with several 
extensions (notably in Secs. PTTDl , PTTTl [P7T1 and [TvT|) . 

We start out in Sec. |J with a discussion of the quantization of the electromagnètic field 
in absorbing or amplifying media. There exists a variety of approaches to this problem 
[9-17]. We will use the method of input-output relations developed by Gruner and Welsch 
|T^,[TB|, and by Loudon and coworkers |T2}- |il||T7j | . The central formula of this section is a 



fluctuation-dissipation relation, that relates the commutator of the operators describing the 
quantum fluctuations in the electromagnètic field to the deviation 1 — S ■ S* from unitarity of 
the scattering matrix S of the system. The relation holds both for absorbing and amplifying 
media. The absorbing medium is in thermal equilibrium at temperature T, and expectation 
vàlues can be computed in terms of the Bose-Einstein function 

f(u,T) = [exp(?iu/k B T)-l]- 1 . (1.1) 

The amplifying medium is not in thermal equilibrium, but the expectation vàlues can be 
obtained from those in the absorbing medium by evaluating the Bose-Einstein function at 
a negative temperature [ P~2"| , |iT • 



In Sec. [IÏ] we apply this general framework to a photodetection measurement. Our 
central result is a relationship between the probability distribution P(n) to count n photons 
in a long time t (long compared to the coherence time of the radiation) and the eigenvalues 
oi, a 2 , . . . cr N of the matrix product S ■ . We call these eigenvalues "scattering strengths". 
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They are between and 1 for an absorbing medium and greater than 1 for an amplifying 
medium. The mean photocount n is proportional to the spectral average N~ x £) n (l — cr n ) of 
the scattering strengths. This spectral average is the absorptivity of the medium, being the 
fraction of the radiation incident in N modes that is absorbed. (An amplifying medium has a 
negative absorptivity) The relation between mean photocount and absorptivity constitutes 
Kirchhoff's law of thermal radiation. We generalize Kirchhoff's law to higher moments of 
the counting distribution by relating the p-th factorial cumulant of n to N" 1 £) n (l — u n ) p - 
While the absorptivity (p = 1) can be obtained from the radiative transfer equation, the 
spectral averages with p > 1 can not. Fortunately, random-matrix theory provides a set of 
powerful tools to compute such spectral averages JO 



We continue in Sec. [TV] with the application of our formula for the photodetection dis- 
tribution to specific random media. We focus on two types of geometries: An open-ended 
waveguide and a cavity containing a small opening. Randomness is introduced by disorder 
or (in the case of the cavity) by an irregular shape of the boundaries. Radiation is emitted 
into iV propagating modes, which we assume to be a large number. (In the case of the cavity, 
N is the number of transverse modes in the opening.) It is unusual, but essential, that all 
the emitted radiation is incident on the photodetector. We show that if only a single mode is 
detected, the counting distribution contains solely information on the absorptivity, while all 
information on higher spectral moments of the scattering strengths is lost. To characterize 
the fluctuations in the photocount we compute the variance Varn = n 2 — n 2 . The variance 
can be directly measured from the auto-correlator of the photocurrent I(t) = I + SI(t), 
according to 

r 00 1 

/ dt6I(0)SI(t) = lim -Varn. (1.2) 

J-oo í^ 00 t 

The bar 7TT indicates an average over many measurements on the same sample. The mean 
photocount n (and hence the mean current I — n/t) contains information on the absorptivity. 
The new information contained in the variance of n (or the auto-correlator of 7) is the 
effective number of degrees of freedom zv eff , defined by |L9| Varn = n(l + n/v e ff)- F° r black- 
body radiation in a narrow frequency interval 5u, one has v e ^ = Nt5u/2ir = v. The counting 
distribution is then a negative-binomial distribution with v degrees of freedom, 

( n _|_ y — l\ 

P{n) oc í 1 exp{-nhu/k B T). (1.3) 

The quantity v e R generalizes the notion of degrees of freedom to radiation from systems that 
are not black bodies. 

A black body has scattering matrix S = 0. (Any incident radiation is fully absorbed.) 
In other words, the scattering strengths a n of a black body are all equal to zero. A ran- 
dom medium, in contrast, has in general a broad (typically bimodal) density of scattering 
strengths. We show that this results in a substantial reduction of below v. In other 
words, the noise in the photocount is anomalously large in a random medium. The reduc- 
tion in í/ e fj holds both for absorbing and amplifying media. The only requirement is a broad 
distribution of scattering strengths. For the random làser, we predict that the ratio v^/u 
vanishes on approaching the làser threshold. No such reduction is expected in a conven- 
tional làser. We discuss the origin of this difference in the concluding Sec. 0, together with 
a discussion of the relationship between z/ eff and the Thouless number of mesoscopic physics. 
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FIG. 1. Scattering geometry consisting of a random medium (dotted) coupled to free space 
via an iV-mode waveguide. The ./V-component vector of outgoing-mode operators a out is linearly 
related to the incoming-mode operators a m and the spontaneous-emission operators b, c. 

II. QUANTIZATION OF THE ELECTROMAGNÈTIC FIELD 

A. Input output relations 

We consider a dielectric medium coupled to free space via a waveguide with N(u) prop- 
agating modes (counting polarizations) at frequency u> (see Fig. [I]). The incoming and out- 
going modes in the waveguide are represented by two iV-component vectors of annihilation 
operators a m (uj), a ont (uj). They satisfy the canonical commutation relations 

K», o£V)] = 5 nm 5(uj - u>), [<», a£(u/)] = 0, (2.1) 
[aT (w), C* V)] = 6 nm 5(u - [a° ut (w), <V)] = 0. (2.2) 

The input -output relations take the form 

a out = 5- a in + C/-6 + y-c + . (2.3) 

The two sets of operators b, and c, é commute with each other and with the set of input 
operators a m ,a in ^. They satisfy the canonical commutation relations 

[b n (u), bl{u')} = 5 nm 5{uü - u'), [b n (u), b m {uü')} = 0, (2.4) 

[ C nM, c m( w ')] = à nm Ó(uJ - üj'), [Cn(u), C m (uj')} = 0, (2.5) 

provided the N x N matrices U(u) and V(u) are related to the scattering matrix S(u) by 

U ■ C/t - v ■ V"< = 1 - S ■ Si (2.6) 

(1 denoting the N x N unit matrix). Equation (|2.6j ) can be understood as a fluctuation- 
dissipation relation: The left-hand side accounts for quantum fluctuations in the electro- 
magnètic field due to spontaneous emission or absorption of photons, the right-hand side 
accounts for dissipation due to absorption (or stimulated emission in the case of an amplify- 
ing medium). Equation ( |2.6|) also represents a link between classical òptics (the scattering 
matrix S) and quantum òptics (the quantum fluctuation matrices U, V). 
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The matrix 1 — S ■ S' is positive definite in an absorbing medium, so we can put V = 
and write 



G out = S ■ a m + U ■ b, U ■ í/ 1 " = 1 - S ■ S ] . (2.7) 

(These are the input-output relations of Ref. |I6|| .) Conversely, in an amplifying medium 
1 — S • is negative definite, so we can put U = and write 

a out = S ■ a m + V ■ c\ V- V* = S ■ - 1. (2.8) 



(The operator c represents the inverted oscillator of Ref. [12|.) Both matrices U,V and 
operators b, c are needed if 1 — 5" • S' is neither positive nor negative definite, which might 
happen if the medium contains both absorbing and amplifying regions. In what follows we 
will not consider that case any further. 



B. Expectation vàlues 

We assume that the absorbing medium is in thermal equilibrium at temperature T. 
Thermal emission is described by the operator b with expectation value 

{bUü>)b m (ü>y = 8 nm 8{u-ü/)f(ü>,T), (2.9) 

where / is the Bose-Einstein function ( |1.1| ). 

The inverted oscillator c accounts for spontaneous emission in an amplifying medium. 
We consider the regime of linear amplification, below the làser threshold. Formally, this 
regime can be described by a thermal distribution at an effective negative temperature —T 



T2| , [ï7j| . For a two-level atòmic system, with level spacing hcuo and an average occupation 
-Wupper > ^V"iower of the two levels, the effective temperature is given by Af upper /iVi ower = 
exp^c^o/feR^)- The zero-temperature limit corresponds to a complete population inversion. 
The expectation value is given by 

(c»4(o/)> = -5 nm 5(cü - u')f{u, — T), (2.10) 



or equivalently by 



c, 



l{uj)c m {uj')) = 5 nm 5{u-u')f{u,T). (2.11) 



(We have used that f{uJ,T) + f(u, —T) = -1.) 

Higher order expectation vàlues are obtained by pairwise averaging, as one would do 
for Gaussian variables, after having brought the operators into normal order (all creation 
operators to the left of the annihilation operators). This procedure is an example of the 
"optical equivalence theorem" |Ï9|,^0l . To do the Gaussian averages it is convenient to 
discretize the frequency as u p = pA, p = 1,2,..., and send A to zero at the end. The 
expectation value of an arbitrary functional JF of the operators b, b^ (or c, c') can then be 
written as a múltiple integral over an array of complex numbers z np [1 < n < N(u p )], 

( : FW>n(fOp)h {blM)} ■■ ) = Z- X J dze-*F[{z np }, {z* np }}, (2.12) 
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with the definitions 



-gè^ (2 - 13) 

Z= fdze-* = Y[ irf{u, >' T) . (2.14) 

J n,p ^ 



The colons in Eq. (|2.12 ) indicate normal ordering, and / dz indicates the integration over 
the real and imaginary parts of all the z np 's. 

III. PHOTODETECTION STATISTICS 
A. General formulas 

We consider the case that the incoming radiation is in the vacuum state, while the 
outgoing radiation is collected by a photodetector. We assume a mode and frequency inde- 
pendent detection efficiency of a photoelectrons per photon. The probability that n photons 
are counted in a time t is given by the Glauber-Kelley-Kleiner formula pï|p2f 

P(n) = -Í(:re- / :>, (3.1) 

I = a f dt'a out Ut')-a out (t'), (3.2) 
Jo 

where we have defined the Fourier transform 

a out (í) = (2tt)- 1/2 / due-' lut a on \üü). (3.3) 
Jo 

The factorial cumulants k v of P(n) are the cumulants of the factorial moments 



n(n — 1) • • • (n — p + 1). For example, k% — n and k 2 = n(n — 1) — n 2 = Varn — n. The 
factorial cumulants have the generating function 

F(0 = E -4- = ln E C 1 + O n P(n) . (3.4) 

p=l P- \n=0 / 

Once is known, the distribution P(n) can be recovered from 

P(n) = — í dzz-^e^-V = hm i-^e F «). (3.5) 

From Eq. ( |3.1|) one finds the expression 

<- r:i: (:<-':). (3.6) 

To evaluate Eq. ( |3.6| ) for the case of an absorbing medium, we combine Eq. ( |2.7| ) with 
Eqs. (|3.2j ) and (|3.3| ), and then compute the expectation value with the help of Eq. ( p.!2j ), 



6 



e F(Ç) = Z 1 J dz exp Z np M np,n' P ' Z n' P ' , 

\ n,p n'.p' ) 

J{UJp,l) Í7T JO 



(3.7) 



(3.8) 



Evaluation of the Gaussian integrals results in the compact expression 

= constant - ln ||M||, (3.9) 

where || • • • || indicates the determinant. (The ^-independent constant can be found from the 
normalization requirement that F(0) = 0.) The matrix U is related to the scattering matrix 
5* by U ■ W = 1 — S ■ [Eq. ( |2.7j )]. This relation determines U up to a transformation 
U — > U ■ A, with A(uj) an arbitrary unitary matrix. Since the determinant ||M|| is invariant 
under this transformation, we can say that knowledge of the scattering matrix sufnces to 
determine the counting distribution. 

The result for an amplifying medium is also given by Eqs. ( |3.8|) and (|3.9| ), with the 



replacement of U by V and f(u p ,T) by —f(uj p ,—T) [in accordance with Eqs. ( |2.8[ ) and 
&))■ 

The determinant ||M|| can be simplified in the limit of large and small counting times 
t. These two regimes will be discussed separately in the next two subsections. A simple 
expression vàlid for all t exists for the mean photocount, 

roc fji ; 

n = t —afTríl-S-S^). (3.10) 
Jo 2tc 

The quantity iV -1 Tr (t — S ■ S<) is the absorptivity, defined as the fraction of the incident 
power that is absorbed at a certain frequency, averaged over all incoming modes. The 
relation ( |3.10| ) between thermal emission and absorption is Kirchhoff's law. It holds also for 
an amplifying medium, upon replacement of f(u,T) by f(u>, — T).Q 



B. Long-time regime 

The long-time regime is reached when uj c t ^> 1, with uj c the frequency interval within 
which S ■ ifr does not vary appreciably. In this regime we may choose the discretization 
Up = pA, A = 2n/t, satisfying 

dí'e^-'V)*' =t5 pp ,. (3.11) 



1 Eq. ([3j]) for an amplifying system is obtained by the replacement of U by V and f(u),T) by 
minus f(uj, —T). Since V ■ V* equals minus 1 — S ■ [Eq. (2.8)], the two minus signs cancel and 
the net result for Eq. ( |3.10|) is that we should replace f(u,T) by plus f(u>, —T). 
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The matrix ( |3.8|) then becomes diagonal in the indices p,p', 



M n P ,ny = s? n ' 5 %s ~ ï a6 PP' ( U ^M ■ U(üj p )) , (3.12) 

t(UJ n , 1) V /nu' 



so that the generating function fl3.9|) takes the form 

F(0 = _í 7o ^Hli-(i-^-^)e«/n 

= ~t ^EM 1 -!!-^»/^?)]. (3.13) 

We have introduced the scattering strengths eri, o"2, . . . erjv, being the eigenvalues of the 
scattering- matrix product S ■ . The result ( |3.13| ) holds also for an amplifying system, 
if we replace f(ui,T) by f(u, —T). 

Expansion of the logarithm in powers of £ yields the factorial cumulants [cf. Eq. ( |3.4j) 1 

Kp =ip-ï)H ^WrEa-^r (3.14) 

Jo 2/r n=1 

The p-th factorial cumulant of P(n) is proportional to the p-th spectral moment of the 
scattering strengths. It is a special property of the long-time regime that the counting 
distribution is determined entirely by the eigenvalues of S ■ S\ independently of the eigen- 
functions. Eq. Q3.14Q can be interpreted as a generalization of Kirchhoff's law (|3.10|) to 
higher moments of the counting distribution. 

C. Short-time regime 

The short-time regime is reached when Q c t <C 1 , with íl c the frequency range over which 
S ■ S* differs appreciably from the unit matrix. (The reciprocal of íl c is the coherence time 
of the thermal emissions.) In this regime we may replace exp[i(cj p — u p i)t'] in Eq. (|3.8|) by 
1, so that M simplifies to 

M np,n' P ' = JT-Tjk ~ ~~o~~ 51 UlmMUmn'M- ( 3 - 15 ) 

J{uJp,l) 2tt m 



If we suppress the mode indices n,n', Eq. (|3.15 ) can be written as 



M p , p , = 77^1 - ^^H) " UM- (3-16) 
The determinant ||M pp /|| can be evaluated with the help of the formulaQ 



2 To verify Eq. ( |3.17| ), take the logarithm of each side and use ln ||M|| = Tr lnM. Then expand 
each logarithm in powers of A and equate term by term. I am indebted to J.M.J. van Leeuwen for 
helping me with this determinant. 
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FIG. 2. A cavity radiates through a hole with an area A. The number N of radiating modes 
at wavelength A is 2irA/\ 2 (counting polarizations). All N modes are detected if the photocathode 
covers the hole. Upon increasing the separation d between hole and photocathode, fewer and fewer 
modes are detected. Finally, single-mode detection is reached when the area of the photocathode 
becomes less than the coherence area ~ d 2 /N of the radiation. 



\\5 pp ,t + A p - B p ,\\ = Hl + ^^-^ll (3.17) 

i 

(with {v4 p }, {B p } two arbitrary sets of matrices). The resulting generating function is[] 

F(0 = -\n\\l-t (1-S-S^af\\. (3.18) 

JO Z7T 

Again, to apply Eq. ( |3.18| ) to an amplifying system we need to replace f(u, T) by f(u>, —T). 

The short-time limit (|3.18|) is more complicated than the long-time limit ( |3.13|) , in the 
sense that the former depends on both the eigenvalues and eigenvectors of S ■ S'. There is 
therefore no direct relation between factorial cumulants of P(n) and spectral moments of 
scattering strengths in the short-time regime. 



D. Single-mode detection 

We have assumed that each of the N radiating modes is detected with equal emciency 
a. At the opposite extreme, we could assume that only a single mode is detected. This 
would apply if the photocathode had an area smaller than the coherence area of the emitted 
radiation (see Fig. Q). Single-mode detection is less informative than multi-mode detection, 
for the following reason. 

Suppose that only a single mode is detected. The counting distribution P(n) is still 
given by Eq. ( |3.1|) , but now / contains only a single element (say, number 1) of the vector 
of operators a out , 

I = a f dt'a°^\t')aTHt'). (3.19) 
Jo 



3 We adopt the convention that the matrix 1 — S ■ S' is embedded in an infinite-dimensional 
matrix by adding zeroes. The matrix 11 outside the integral over ui in Eq. ( 3.18| ) is then interpreted 
as an infinite-dimensional unit matrix. 
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L 

FIG. 3. Disordered waveguide (length L) connected at one end to a photodetector. 



This amounts to the replacement of the matrix U in Eq. ( |3.8D by the projection V ■ U, with 
V = 5 nm 5 n i. Instead of Eq. ( |3.10| ) we have the mean photocount 



n 



, dndntaf r .^±,, , 



, 00 duo ( 2n£ dn\ , 
F(0 = -t 7^ In 1— (3.21) 



The generating function now takes the form 

r°°d(jj , 
/ — h- 

lo 2n \ t duj 
in the long-time regime, and 

F(0 = -ln(l-£n) (3.22) 

in the short-time regime. We see that the entire counting distribution is determined by 
the mean photocount, hence by the absorptivity of the single detected mode. This strong 



version of Kirchhoff's law is due to Bekenstein and Schiffer f23|| . It holds only for the case 
of single-mode detection. Multi- mode detection is determined not just by n, being the flrst 
spectral moment of the scattering strengths, but also by higher moments. 



E. Waveguide geometry 

Figs. [I] and show a cavity geometry. Alternatively, one can consider the waveguide 
geometry of Fig. || The waveguide has cross-section A, corresponding to N = 2-nA/X 2 
modes at wavelength A. The 2N x 2N scattering matrix S consists of four N x N blocks, 

S =([,;), (3-23) 

namely two reflection matrices r, r' (reflection from the left and from the right) and two 
transmission matrices t, t' (transmission from right to left and from left to right). Reciprocity 
relates t and t' (they are each others transpose). 

A photodetector detects the radiation emitted at one end of the waveguide, while the 
radiation emitted at the other end remains undetected. If the radiation from both ends 
would be detected (by two photodetectors), then the eigenvalues of S ■ would determine 
the counting distribution in the long-time limit, as in the cavity geometry. But for detection 
at the left end only, one needs instead the eigenvalues of the matrix r-r^+í-í^ (or r'-r'' +t'-t'> 
for detection at the right end). More precisely, the general expression for the characteristic 
function is given by Eqs. ( |3.7| ) and ( |3.8| ) upon replacement of the 2N x 2N matrix U by the 
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projection V ■ U, with V nm = lifl<n = m<iV and V nm = otherwise. In the long-time 
regime one obtains 

fOO Al ] 

F(0 = ~t —ln\\l-(l-r-rt-t-t^af\\, (3.24) 

JO 27T 

and in the short-time regime 

F{£) = -ln||l -t / — (1 -r ■ r f -t-tf)Çaf\\. (3.25) 
Jo 2n 

The eigenvalues of r-r'+t-t' differ from the sum R n +T n of the reflection and transmission 
eigenvalues (eigenvalues oïr • and t ■ t\ respectively), because the two matrices r • and 
t ■ t ' do not commute. The absorptivity 

N 

N~ l Tr (t-r -r ] -t-t ] ) = N- 1 ^^- Rn~ T n ) 

n=l 

does depend only on the reflection and transmission eigenvalues. It determines the mean 
photocount 



roc Ai ; 

t — a f Tr (U - r • r f - í ■ t f ), (3.26) 

JO 27T 



n 



in accordance with Kirchhoff's law. Higher moments of the counting distribution can not 
be obtained from the reflection and transmission eigenvalues, but require knowledge of the 
eigenvalues of r ■ + t ■ t\ A substantial simplification occurs, in the case of absorption, if 
the waveguide is sufnciently long that there is no transmission through it. Then t ■ t* can be 
neglected and the counting distribution depends entirely on the reflection eigenvalues. 



IV. APPLICATIONS 
A. Black-body radiation 



Let us first check that we recover the familiar results for black-body radiation ||l9| , |20 
The simplest "step-function model" of a black body has 





for 


— Lüç)\ 


c 


for 


\üJ — üüo\ 



Incident radiation is fully absorbed within the frequency interval fl c around c^o and fully 
reflected outside this interval. Typically, Q c ojq so that we may neglect the frequency 
dependence of N(uS) and f(u,T), replacing these quantities by their vàlues at uj — uj q . 
The generating function ( |3.13| ) in the long-time regime then becomes 

F{t) = -^Hl-t"f)- (4-2) 
The inversion formula (|3.5[) yields the counting distribution 
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1 ' n\T(v) (l + n/vY+v- K ' 

This is the negative-binomial distribution with v = Ntfl c /2^ degrees of freedom. [For 
integer v, the ratio of Gamma functions forms the binomial coeficient that counts 

the number of partitions of n bosons among v states, cf. Eq. (|1.3|) .] Note that v 1 
in the long-time regime. The mean photocount is n = vaf. In the limit n/v — > 0, the 
negative-binomial distribution tends to the Poisson distribution 

P(n) = \n n e~ n (4.4) 
ni 

of independent photocounts. The negative-binomial distribution describes photocounts that 
occur in "bunches" . Its variance 

Var n = n(l + n/v) (4-5) 

is larger than the Poisson value by a factor 1 + n/v. 
Similarly, the short-time limit ( [3.18|) becomes 



F(£) = -iVln (l - ^«/) , (4-6) 

corresponding to a negative-binomial distribution with iV degrees of freedom. 

In the step-function model (|4.1|) S changes abruptly from to 1 at \uj — ujq\ = |Í2 C . 
A more realistic model would have a gradual transition. A Lorentzian frequency profile is 
commonly used in the literature |p4}| , for the case of single-mode detection. Substitution of 
1 - (S ■ St)n = [1 + A(lü - luq) 2 /^}- 1 into Eq. (^2Üj) and integration over u in Eq. flgT2l 



(neglecting the frequency dependence of /) yields the generating function in the long-time 
regime, 



F(£) = kíl c ll-Jl-taf). (4.7) 



The corresponding counting distribution is 



Jl + af) , (4.8) 



with n = \tíl c af and K a Bessel function. [The normalization constant is C = 
exp(|tí2 c )(tí7 c /7r) 1 / 2 (l + af) 1 ^.] This distribution was first obtained by Glauber pï| . It 
is closely related to the socalled ií-distribution in the theory of scattering from turbulent 
media p5| , ^6| . The counting distribution ( |3.22|) in the short-time regime remains negative- 



binomial. 

In most realizations of black-body radiation the value of the Bose-Einstein function 
f(uo,T) is <^ 1. The difference between the two distributions ( |4.3| ) and ( |4.8| ) is then quite 



small, both being close to the Poisson distribution ( |4.4|) . 
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B. Reduction of degrees of freedom 



We now turn to applications of our general formulas to specific random media. We 
concentrate on the long-time regime and assume a frequency-resolved measurement, in which 
photons are only counted within a frequency interval 5uj around c^o- (For a black body, this 
would correspond to the step-function model with Q c replaced by Su.) We take óu smaller 
than any of the characteristic freqüències u c , Q c , but necessarily greater than 1/t. The 
factorial cumulants are then given by 

1 N 

K p =(p-l)\v(afy-Y,a-VnY, (4.9) 

iV n=l 

where v = Nt5u/2Tí. For comparison with black-body radiation we parameterize the vari- 
ance in terms of the effective number of degrees of freedom |T9[ , 



Varra = n(l +n/v eS ), (4.10) 
with u e g = v for a black body [cf. Eq. (|4.5|) ]. Eq. (|4.9| ) implies 

^eff [Ení 1 - <7n)Y 



V iVEn(l-^n) 2 



< I- (4.11) 



We conclude that the súper- Poissonian noise of a random medium corresponds to a black 
body with a reduced number of degrees of freedom. The reduction occurs only for multi- 
mode emission. (Eq. ( [4 . 1 1| ) with N — 1 gives z/ e g = v.) In addition, it requires multi- 



mode detection to observe the reduction, because single-mode detection contains no other 
information than the absorptivity (cf. Sec. [111 D| ). 

An ensemble of random media has a certain scattering-strength density 



p{a) = ly j 8{a-a n )), (4.12) 




where the brackets (■ ■ •) denote the ensemble average. In the large-iV regime sample-to- 
sample fluctuations are small, so the ensemble average is representative for a single system.0 
We may therefore replace J2 n by / da p(cr) in Eqs. (|4.9| ) and fl4.11|) . In the applications 



that follow we will restrict ourselves to the large-iV regime, so that we can ignore sample- 
to-sample fluctuations. All that we need in this case is the function p(a). Random-matrix 
theory |TS| provides a method to compute this function for a variety of random media. 



4 This statement is strictly speaking not correct for amplifying systems. The reason is that the 
ensemble average is dominated by a small fraction of members of the ensemble that are above 
the làser threshold, and this fraction is non-zero for any non-zero amplification rate. This is a 
non-perturbative finite-iV effect that does not appear if the ensemble average is computed using 
the large-A^ perturbation theory employed here. 
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C. Disordered waveguide 



As a first example we consider the thermal radiation from a disordered absorbing wave- 
guide (Fig. |J). The length of the waveguide is L, the transport mean free path in the medium 
is l, the velocity of light is c, and r s = l/c is the scattering time. The absorption time r a is 
related to the imaginary part e" > of the (relative) dielectric constant by l/r a = uj e". We 
assume that r s and r a are both l/u , so that scattering as well as absorption occur on 
length scales large compared to the wavelength. We define the normalized absorption ratef] 

7 = ^- (4.13) 

We call the system weakly absorbing if 7 <C 1, meaning that the absorption is weak on the 
scale of the mean free path. If 7 ^> 1 we call the system strongly absorbing, 7 — > 00 being 
the black-body limit. For simplicity we restrict ourselves to the case of an infinitely long 
waveguide (more precisely, L ^> so that transmission through it can be neg lected.Q 

In the absence of transmission the scattering matrix S coincides with the reflection matrix 
r, and the scattering strengths a n coincide with the reflection eigenvalues R n (eigenvalues 
of r • r>). The density p(a) for this system is known for any value of N [27,p8||. The general 



expression is a series of Laguerre polynomials, which in the large- N regime of present interest 
simplifies to 

m = < „ < * . (4,4) 

7T 1 + 17 

(The large- N regime requires N ^> 1, but in weakly absorbing systems the condition is 
stronger: iV ^> I/-/7.) This leads to the effective number of degrees of freedom 

V = ^^S = ^ + ^ + d + V7)-/V. (4.15) 

plotted in Fig. [|, with a mean photocount of 

n = \uah (v'l + 4/7 - l) • (4.16) 

For strong absorption, 7 > 1, we recover the black-body result u cS = u, as expected. For 
weak absorption, 7 <C 1, we find z/ efï = 2vy/j. 



5 The coeficient 16/3 in Eq. ( 4.1 3| ) is chosen to facilitate the comparison between waveguide 
and cavity in the next subsection, and refers to three-dimensional scattering. In the case of two- 
dimensional scattering the coefficient is ir 2 /2. The present definition of 7 differs from that used in 



Ref. [27| by a factor of two. 



6 Results for an absorbing waveguide of finite length follow from Eqs. ( 4.3C| )-( [4.32| ) upon changing 
the sign of 7. 
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black-body limit 




amplification 

làser threshold 



FIG. 4. Effective number of degrees of freedom as a function of normalized absorption or 
amplification rate. The dashed curve is Eq. ( 4. 15j ) for an absorbing, infinitely long disordered 
waveguide, the sòlid curves are Eqs. ( 4,24 ) and ( |4.29| ) for the chaotic cavity. For the cavity both 
the cases of absorption and amplification are shown. (See Fig. ^ for the amplifying waveguide.) 
The black-body limit for absorbing systems and the làser threshold for amplifying systems are 
indicated by arrows. 



The characteristic function in the long-time frequency-resolved regime follows from 

F(0 = ~jjj dcrp(a) ln[l - (1 - a)Çaf}. (4.17) 

Substitution of Eq. ( |4.14j ) into Eq. fl4.17|) yields a hypergeometric function, which in the 
limit 7 <C 1 of weak absorption simplifies to 

F(0 = u eS (l-y/l- &f) , v BS = 2uy/ï. (4.18) 
The counting distribution corresponding to Eq. (|4.18|) , 

P(n) oc i í^=j) K n-i/2 (vesy/ï+^f), (4-19) 

is Glauber's distribution ( |4.8j ) with an effective number of degrees of freedom. Note that Eq. 
( |4.8|) resulted from single-mode detection over a broad frequency range, whereas Eq. (|4.19|) 
results from multi-mode detection over a narrow frequency range. It appears as a coincidence 
that the two distributions have the same functional form (with different parameters). 



D. Chaotic cavity 

Our second example is an optical cavity radiating through a small hole covered by a 
photodetector (Fig. H). The area A of the hole should be small compared to the surface 
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0.5 1 

scattering strength a 

FIG. 5. Sòlid curves: Scattering-strength density of an absorbing chaotic cavity in the large- iV 
regime, calculated from Eq. (5.6) for four vàlues of the dimensionless absorption rate 7. The density 



( 4.14 ) for an absorbing, infinitely long disordered waveguide is included for comparison (dashed). 
For 7>1 the results for cavity and waveguide coincide. 



area of the cavity. The cavity should have an irregular shape, or it should contain random 
scatterers — to ensure chaotic scattering of the radiation inside the cavity. It should be 
large enough that the spacing Au of the cavity modes near frequency u is <C uq. For this 
system we define the normalized absorption rate as 



Tdwell _ _ 2n 
7 = ; r dwell = 



NAu 



(4.20) 



The time Td we n is the mean dwell time of a photon in the cavity without absorption. The 
frequency l/r^n represents the broadening of the cavity modes due to the coupling to 
the N = 2-iïA/X 2 modes propagating through the hole. The broadening is much greater 
than the spacing Au for N ^> 1. The large- N regime requires in addition iV ^> I/7. The 
scattering-strength density in the large- N regime can be calculated using the perturbation 
theory of Ref. ^ . 



The result is a rather complicated algebraic function, see the Appendix. It has a simple 
form in the limit 7 <C 1 of weak absorption, 



N 



(7 



,1/2 



fa 



1/2 



2tt 



o~_ < a < a + , 



(4.21) 



with a± = 1 — 37 ± 27V / 2· In the opposite limit 7 ^> 1 of strong absorption, p(a) is given 
by the same Eq. Q4.14j ) as for the infinitely long disordered waveguide. The crossover from 
weak to strong absorption is shown in Fig. [5]. The value 7 = 1 is special in the sense that 
p(a) goes to zero or infmity as o — > 0, depending on whether 7 is smaller or greater than 1. 
We find the mean and variance of the photocount 
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va f 7 . 

n = —11, 4.22 
1 + 7 

Varn = ra + i/( a /)V ^7r-^T- , ( 4 -23) 



;i+7) 

corresponding to the effective number of degrees of freedom 

= (I + 7) 2 
i/ 7 2 + 27 + 2 ' 



(4.24) 



Again, v c ff = v for 7 ^> 1. For 7 < 1 we find f e g = This factor-of-two reduction of the 
number of degrees of freedom is a "universal" result, independent of any parameters of the 
system. The chaotic cavity is compared with the disordered waveguide in Fig. [|. The ratio 
v c <s_jv for the chaotic cavity remains finite no matter how weak the absorption, while this 
ratio goes to zero when 7 — > in the case of the infinitely long disordered waveg uide.Q 

E. Random làser 

The examples of the previous subsections concern thermal emission from absorbing sys- 
tems. As we discussed in Sec. |111 BL our general formulas can also be applied to amplified 



spontaneous emission, by evaluating the Bose-Einstein function / at a negative tempera- 



ture [p]2])[17| . Complete population inversion corresponds to / = — 1. The amplification rate 
l/r a = c^ok"| should be so small that we are well below the làser threshold, in order to 
stay in the regime of linear amplification. The làser threshold occurs when the normalized 
amplification rate 7 reaches a critical value j c . (Sample-to-sample fluctuations in the làser 
threshold |30| are small in the large-iV regime.) For the cavity 7 C = 1. For the disordered 



waveguide one has 

7c = {jpj if L » l. (4.25) 

Since 7 C — > in the limit L — > 00, the infinitely long waveguide is above the làser threshold 
no matter how weak the amplification. 

A duality relation [^] between absorbing and amplifying systems greatly simplifies the 
calculation of the scattering strengths. Dual systems differ only in the sign of the imaginary 
part e" of the dielectric constant (positive for the absorbing system, negative for the am- 
plifying system). Therefore, dual systems have the same value of r a and 7. The scattering 
matrices of dual systems are related by Si = hence SL • SL. = (S + ■ S+)~L (The 
subscript + denotes the absorbing system, the subscript — the dual amplifying system.) 
We conclude that the scattering strengths ai, 02, . . . cxtv of an amplifying system are the 
reciprocal of those of the dual absorbing system. The densities p±(cr) are related by 



7 For a waveguide of finite length L (S> l) one has instead v e s/v — > 5//L in the limit 7^0 (cf. 
Sec. pVH). 



17 




1 1.5 2 

scattering strength a 

FIG. 6. Scattering-strength density of an amplifying chaotic cavity, calculated from the results 
in Fig. H by means of the duality relation ( 4.26| ). First and higher moments diverge for 7 > 1. 



p + (l/a). 



(4.26) 



In Fig. ^ we show the result of the application of the transformation ( |4.26|) to the densities 
of Fig. [| for the case of a cavity. The critical value 7 C = 1 is such that the first and higher 
moments are finite for 7 < 7 C and infinite for 7 > 7 C . 

We find that the expressions for n, V&rn, and u e g/u in the amplifying chaotic cavity 
differ from those in the dual absorbing cavity by the substitution of 7 by —7: 



n 



vafj_ 

'i-r 



7 

Varn = n + u{af) 2 ^ 2 — 



2 7 + 2 



1Y 



V 



(1-7) 1 



r 



27 + 2 



(4.27) 
(4.28) 
(4.29) 



The Bose-Einstein function / is now to be evaluated at a negative temperature, so that / < 0. 
In Fig. § we compare v e s / v fo f amplifying and absorbing cavities. In the limit 7 — > the two 
results coincide, but the 7-dependence is strikingly different: While the ratio v e g/v increases 
with 7 in the case of absorption, it decreases in the case of amplification — vanishing at 
the làser threshold. Of course, close to the làser threshold [when 7 > 1 — (fi c Tdwcii) -1 ^ 2 ] the 
approximation of a linear amplifier breaks down and a non-linear treatment (along the lines 
of Ref. P2"| ) is required. 

For the amplifying disordered waveguide we can not use the infinite-length formulas of 
Sec. [IV (J| , because then we would be above threshold for arbitrarily small 7. The counting 
distribution P(n) at a finite L requires the density of eigenvalues of the matrix r • + t ■ t>, 
as explained in Sec. |III E| . The density itself is not known, but its first few moments have 
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0.2 0.4 0.6 0.8 1 

l/lc 

FIG. 7. Effective number of degrees of freedom for an amplifying disordered waveguide of 
finite length L (much greater than the transport mean free path l), computed from Eq. ( [4.32 ). 
The làser threshold occurs at 7 = j c = (47rZ/3L) 2 . 



been calculated recently by Brouwer [ 33fl . This is suficient to compute z/ c fj, since we only 
need the first two moments of P(n). The results for 7, 7 C < 1 areQ 



n 



sm s 



cos s 



Vam = n + v(af) 



V 



2 sin 4 s 
— 3 sin 3 s — cos 3 s sin s 

2^7(1 



s — s cos s + s sin s + sin s 



coss) 2 sin 2 s 



s — s cos s + s sin s + sin s — 3 sin s — cos 3 s sin s 



(4.30) 

(4.31) 
(4.32) 



where we have abbreviated s = ny'j/'jc. Fig. |7| shows a plot of Eq. (|4.32|) . Notice the 
limit v e ïï/v = 51/ L for r y/ r y c — > 0. The reduction of the number of degrees of freedom on 
approaching the làser threshold is qualitatively similar to that shown in Fig. ^ for the chaotic 
cavity. 



F. Broad-band detection 



In these applications we have assumed that only photons within a narrow frequency 
interval óu are detected. This simplifies the calculations because the frequency dependence 



8 Ref. considers an absorbing waveguide. The amplifying case follows by changing the sign of 
the parameter 7. Eq. (13c) in Ref. [33 contains a misprint: The second and third term between 
brackets should have, respectively, signs minus and plus instead of plus and minus. 
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of the scattering matrix need not be taken into account. In this subsection we consider the 
opposite extreme that all freqüències are detected. We will see that this case of broad-band 
detection is qualitatively similar to the case of narrow-band detection considered so far. 
We take a Lorentzian frequency dependence of the absorption or amplification rate, 



7(w) 



7o 



i + 4(cu-^ ) 2 /r 2 ' 

The characteristic frequency Q c for the scattering strengths is defined by 



iyi+70- 



(4.33) 



(4.34) 



The two freqüències íl c and T are essentially the same for 70 ~ 1, but for 70 ^> 1 the former 
is much bigger than the latter. The reason is that what matters for the deviation of the 
scattering strengths from zero is the relative magnitude of 7(0;) with respect to 1, not with 
respect to 70 . As in Sec. |V A| , we assume that Q c <C u , so that we may neglect the frequency 
dependence of N and /. The mean and variance of the photocount in the long-time regime 
are given by 



n 



taf 



du 
2^ 



da p(a, u)(l 
du 



a 



— / da p(a, üj)(1 
2tt J 



- cr 



(4.35) 
(4.36) 



Again, we have assumed that N is sumciently large that sample-to-sample fluctuations can 
be neglected and we may replace J2 n by J da. The scattering-strength density p depends on 
ui through the rate 7(0;). 

In the absorbing, infinitely long disordered waveguide p oc w^y for 7 <C 1, hence the 
integrands in Eqs. ( |4.35| ) and (|4.36|) decay oc 1/\uj — lúq\ and the integrals over u diverge. 
A cutoff is provided by the finite length L of the waveguide. When 7 drops below (l/L) 2 , 
radiation can be transmitted through the waveguide with little absorption. Only the fre- 
quency range \u — uq\ £ T(L/l)^yo, therefore, contributes effectively to the integrals ([4.35 



and (^4.361) . To leading order in (L/í)^/7o we can take the infinite-L result for p with the 
cutoff we mentioned. The result (for {l/L) 2 <C 70 1) is 



NtT 



n 



2tt 



ln 



Var n 



L 
1 

NtT , 



0(1) 



ln( T V^)+C(l) 



If we write Varn = n(l + n/u e g), as before, then 



NtT 



7T 



ln(yV7^)+C(l) 



(4.37) 
(4.38) 

(4.39) 



In the case of narrow-band detection considered in Sec. [IV i] we had u^/Ntòuj = ix "^70 
for (l/L) 2 <C 7o 1. The difference with Eq. ( |4.39| ) (apart from the replacement of 8uj by 
T Pà Q c ) is the logarithmic enhancement factor, but still v e g NtT for 70 <C 1. 



20 



0.6 



0.4 



0J 



0.2 




FIG. 8. Plot of Eq. (4.42) for the effective number of degrees of freedom of an absorbing 
or amplifying chaotic cavity, in the case of broad-band detection with a Lorentzian frequency 
dependence 7(0») = 70 [1 + 4(w — uq) 2 /T 2 ]~ l of the absorption or amplification rate. For 7q < 1 the 
curves are qualitatively similar those plotted in Fig. || for the case of narrow-band detection. In 
the absorbing cavity v e s increases oc -y/73 as 70 — > 00, instead of saturating as in Fig. ||, because 
the characteristic frequency fï c increases oc in that limit [cf. Eq. ( [4.34 )]. 



For the chaotic cavity, we can compute n and Var n directly from the narrow-band results 
Q ), (Ül), ÇgT27|) , and (|428| ), by substituting Eq. (g3g ) for 7 and integrating over u. 



There are no convergence problems in this case. The results are 

fi = ±Aríra/ vfe' (440) 

V^n + Nnaf) ^ 2 ^ , (4.41) 

" - 4(1 ± 70)6/2 (4.42) 



NtT 9 7 ^ ± 2O70 + 16 

The ± indicates that the plus sign should be taken for absorption and the minus sign for 
amplification. The function ( |4.42| ) is plotted in Fig. In the strongly absorbing limit 



7o — > 00, the effective number of degrees of freedom z/ e g — > |iVíf2 c , which up to a numerical 
coefficient corresponds to the narrow-band limit z/ e g — > Nt5uj/2n upon replacement of 8uj 
by Q c . The limit 70 — > is the same for absorption and amplification, z/ efT — > ^Ntíl c , again 
corresponding to the narrow-band result z/ e ff — > NtSu/Air up to a numerical coefficient. 
Finally, the ratio Ueg/NtT tends to zero upon approaching the làser threshold 70 — > 1 in an 
amplifying system, similarly to the narrow-band case. The qualitative behavior of v^jv is 
therefore the same for broad-band and narrow-band detection. 
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V. CONCLUSION 



A. Summary 

In conclusion, we have shown that the photodetection statistics contains substantially 
more information on the scattering properties of a medium than its absorptivity. The 
mean photocount n is determined just by the absorptivity as dictated by Kirchhoff's law. 
Higher order moments of the counting distribution, however, contain information on higher 
spectral moments of the scattering strengths a n (being the N eigenvalues of the scattering 
matrix product S • S^). These higher moments are independent of the absorptivity, which is 
determined by the first moment. While the absorptivity follows from the radiative transfer 
equation, higher spectral moments are outside of the range of that approach. We have used 
random-matrix theory for their evaluation. 

To measure these higher spectral moments, it is necessary that the counting time t 
is greater than the coherence time l/íí c of the thermal radiation (being the inverse of the 
absorption or amplification line width). It is also necessary that the area of the photocathode 
is greater than the coherence area (being the area corresponding to one mode emitted by the 
medium). Single-mode detection yields solely information on the absorptivity. Multi-mode 
detection is unusual in photodetection experiments, but required if one wants to go beyond 
Kirchhoff's law. 

We have shown that the variance Varn of the photocount contains information on the 
width of the density p(er) of scattering strengths. We have computed this density for a 
disordered waveguide and for a chaotic cavity, and find that it is very wide and strongly non- 
Gaussian. In an absorbing medium, the deviations from Poisson statistics of independent 
photocounts are small because the Bose-Einstein function is 1 for all practical freqüències 
and temperatures. Since the Poisson distribution contains the mean photocount as the only 
parameter, one needs to be able to measure the super-Poissonian fluctuations in order to 
obtain information beyond the absorptivity. The deviations from Poisson statistics are easier 
to detect in an amplifying medium, where the role of the Bose-Einstein function is played 
by the relative population inversion of the atòmic states. 

We have shown that the super-Poissonian fluctuations in a linearly amplifying random 
medium are much greater than would be expected from the mean photocount. If we write 
Varn — n = n 2 /V e ff, then v c s would equal Nt5uj/2n = v if all N modes reaching the 
photodetector would have the same scattering strength. (We assume for simplicity that 
only a narrow band óu is detected, in a time í; For broad-band detection 5cu should be 
replaced by Q c .) The effective number z/ c g of degrees of freedom is much smaller than 
v for a broad p(cr), hence the anomalously large fluctuations. On approaching the làser 
threshold, the ratio v e g / v goes to zero. In a conventional làser the noise itself increases with 
increasing amplification rate because n increases, but v e s does not change below the làser 
threshold. Typically, light is emitted in a single mode, hence u e g equals t5uj/2-n independent 
of the amplification rate. In a random làser a large number N of cavity modes contribute 
to the radiation, no matter how small the frequency window Sou, because the cavity modes 
overlap. The overlap is the consequence of the much weaker confinement created by disorder 
in comparison to that created by a mirror. The reduction of the number of degrees of 
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freedom is a quantum optical effect of overlapping cavity modes that should be observable 
experimentally. 



B. Relation to Thouless number 



The Thouless number iVr plays a central role in mesoscopic physics [34 . It is a dimen 



sionless measure of the coupling strength of a closed system to the outside world, 

r dwcllAu; ^ 

(We use ~ instead of = because we are ignoring numerical coefficients of order unity.) As 
before, Au; is the spacing of the eigenfrequencies of the closed system and r dwell is the mean 
time a particle (electron or photon) entering the system stays inside. In a conducting metal, 
ÍVt is the conductance in units of the conductance quantum e 2 /h. The metal-insulator 
transition occurs when becomes of order unity. It is assumed that there is no absorption 
or amplification, as is appropriate for electrons. 

For the two types of systems considered in this work, one has ÍVt — Nl/L for the 
disordered waveguide and ~ iV for the chaotic cavity. We notice that is related to 
the effective number of degrees of freedom in the limit of zero absorption and amplification, 

hm ^ (5.2) 

The ratio of v e g to the black-body value v is the same as that of A't to the number of 
propagating modes N. We believe that the relation (|5.2| ) holds for all random media, not 
just for those considered here. 
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APPENDIX. SCATTERING-STRENGTH DENSITY OF A CHAOTIC CAVITY 

As derived in Ref. ||35|1 , absorption in a chaotic cavity (with rate l/r a ) is statistically 
equivalent to the loss induced by a fictitious waveguide that is weakly coupled to the cavity. 
The coupling has transmission probability T' for each of the N' modes in the fictitious wave- 
guide. The equivalence requires the limit N' — > oo, V — > 0, at fixed N'V = 27r/r a Au; (with 
Au; the spacing of the cavity modes). 

We are therefore led to consider the system illustrated in Fig. BJ A chaotic cavity without 
absorption containing two openings. One opening is coupled to an A-mode waveguide with 
transmission probability 1 per mode, the other opening is coupled to the fictitious A'-mode 
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FIG. 9. An absorbing cavity with one opening is statistically equivalent to the non-absorbing 
cavity with two openings shown here. The thin line in the second opening indicates a barrier with 
transmission probability V for each of the N' modes in the waveguide attached to the opening. 
(The limit r' — > 0, N' — > oo at fixed N'T' is required for the equivalence with absorption.) 



waveguide with transmission probability T' per mode. The scattering strength a n equals 
1 — T n , with T n an eigenvalue of the transmission-matrix product t ■ v (t being the N x N' 
transmission matrix from one waveguide to the other). The density of transmission eigenval- 
ues p(T) = (S n 5(T — T n )) can be calculated in the large-iV regime using the perturbation 
theory of Ref. p9[ . The scattering-strength density p(a) then follows from a = 1 — T. 
The result is non-zero for <r m i n < a < a max , with the definitions 



cr_ if 7 < 1, 
if 7 > 1, 



(Tu 



<y± 



8 + 20 7 2 -7 4 ± 7 (8 + 7 2 ) 3 / 2 
8(l + 7) 3 ' 



(5.3) 
(5.4) 
(5.5) 



[We use the same dimensionless absorption rate 7 = NT'/N = 2tt/Nt 3ü Auj as in Eq. (|4.20 ). 
Inside this interval the scattering-strength density is given by 



7T 



a + b) 



1/3 



1 2 



27 + 2 - 6a 



i-b) 1 '* 
[a + b) 1 / 3 - 



[a + b) 1 / 3 
(a-6) 1 / 3 



+ (a - 

2\ -1 



b) 



1/3 



a 
b 



( 7 -l) 3 + 9(l + | 7 2 )^ 



(3 + 3 7 ) 3/2 [a(a- ( x_)(a + -a)] 1 / 2 . 
Eq. ( |5.6|) is plotted in Fig. |5] for several vàlues of 7. 



(5.6) 

(5.7) 
(5.8) 
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